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We have searched for gravitational waves from coalescing low mass compact binary systems with a total 
mass between 2 and 35 Mq and a minimum component mass of 1 Mq using data from the first year of the 
fifth science run (S5) of the three LIGO detectors, operating at design sensitivity. Depending on mass, we are 
sensitive to coalescences as far as 150 Mpc from the Earth. No gravitational wave signals were observed above 
the expected background. Assuming a compact binary objects population with a Gaussian mass distribution 
representing binary neutron star systems, black hole-neutron star binary systems, and binary black hole systems, 
we calculate the 90%-confidence upper limit on the rate of coalescences to be 3.9 x 10 -2 yr -1 !/^ 1 , 1.1 x 
10 -2 yr _1 L^ 1 , and 2.5 x 10~ 3 yr~ 1 L ] ^ 1 respectively, where Lio is 10 10 times the blue solar luminosity. We 
also set improved upper limits on the rate of compact binary coalescences per unit blue-light luminosity, as a 
function of mass. 

PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.60.Jd, 97.60.Lf, 97.80.-d 



I. INTRODUCTION 

Among the most promising candidates for the first detec- 
tion of gravitational waves (GW) are signals from compact bi- 
nary coalescences (CBC), which include binary neutron stars 
(BNS), binary black holes (BBH), and black hole-neutron star 
binaries (BHNS). The inspiral waveforms generated by these 
systems can be reliably predicted using post-Newtonian (PN) 
perturbation theory, until the last fraction of a second prior to 
merger. These waveforms can be used in matched filtering of 
noisy data from gravitational wave detectors to identify GW 
candidate events. 

Astrophysical estimates for CBC rates depend on a number 
of assumptions and unknown model parameters, and are still 
uncertain at present. In the simplest models, the coalescence 
rates should be proportional to the stellar birth rate in nearby 
spiral galaxies, which can be estimated from their blue lumi- 
nosity [ 1]; we therefore express the coalescence rates per unit 
Liq, where Lio is 10 10 times the blue solar luminosity (the 
Milky Way contains - 1.7L 10 |Hl). The most confident BNS 
rate predictions are based on extrapolations from observed bi- 
nary pulsars in our Galaxy; these yield realistic BNS rates of 
5 x 1CP 5 yr^L^ 1 , although rates could plausibly be as high 
as 5 x 1CT 4 yr" 1 !,^ 1 fill- Predictions for BBH and BHNS 
rates are based on population synthesis models constrained 
by these and other observations. Realistic rate estimates are 
2 x 10" 6 yr" 1 ^ 1 for BHNS H and 4 x 1CT 7 vr^L^ 1 for 
BBH U; both BHNS and BBH rates could plausibly be as 
high as 6 x 1(T 5 yr- 1 /^ 1 g@|. 

The Laser Interferometer Gravitational-wave Observatory 
(LIGO) detectors achieved design sensitivity in 2005, and 
completed a two-year-long science run (S5) in November 
2007. Results from searches for GW from CBC by the 



LIGO Scientific Collaboration (LSC) using data from previ- 
ous science runs with ever-increasing sensitivity are reported 
inRefs. SB HUSH! ■ 

This paper summarizes the search for GW signals from 
CBC with component masses greater than or equal to 1 so- 
lar mass (Mq) and total mass ranging from 2 to 35 Mq, using 
the first year of data from LIGO's S5 run, between Novem- 
ber 4th, 2005 and November 14th, 2006. During this time, 
the LIGO detectors were sensitive to signals from CBC with 
horizon distances (Table HB of 30 Mpc for BNS (25 seconds 
in the LIGO detection band) and 150 Mpc for systems with a 
total mass of ~ 28 Mq (0.5 seconds in the LIGO band). Sub- 
sequent papers will report the results of similar searches using 
the data from the second year (during which time the Virgo 
detector was in observational mode), searches for higher-mass 
systems (between 25 and 100 Mq), and specialized searches 
targeting particular subsets of signals. 

The component objects of true astrophysical compact bina- 
ries will in general have some angular momentum, for which 
PN waveforms that incorporate non-zero values for the spin 
parameters are available 11121 fl3ll . However, for most of the 
parameter space, the effect of spin on the waveforms is small, 
and the signals can be captured using non-spinning waveform 
templates (Appendix [1} with only a small loss in the signal- 
to-noise ratio (SNR); this is the approach taken in the search 
described here. In some regions of parameter space, the effect 
of spin is larger, and dedicated searches lfl3l [l4l [l5l [l6Tl may 
be more effective. The LSC continues to develop more effec- 
tive methods for searching for signals with strong modulations 
due to spin. 

The rest of this paper is organized as follows. Section HT1 
summarizes the search pipeline that was employed. Section 
Hill describes the output of the search: detection candidate 
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events which are examined and rejected using a detection con- 
fidence procedure. Section|IV]describes the evaluation of our 
detection efficiency using simulated GW signals injected into 
the detectors' data streams. Section [V] discusses the upper 
limit calculation that was performed, and the resulting upper 
limits on CBC rates neglecting spin of the coalescing objects. 
SectionlVTldiscusses how our sensitivity is affected when spin 
is included. Finally, Section lVnl presents the conclusions, fol- 
lowed by several appendices on certain technical aspects of 
the search. 



II. THE DATA ANALYSIS PIPELINE 

The pipeline used for this analysis has been described in 
previous documents 

OH dim 

and was used to search 
for BNS in LIGO's third and fourth science runs Hill . The 
main aspects of the pipeline and new features used here are 
detailed below. 

The data analysis proceeds as follows. The gravita- 
tional wave strain data are recorded from each of the three 
LIGO detectors: the HI and H2 detectors at LIGO Hanford 
Observatory (LHO) and the LI detector at LIGO Livingston 
Observatory (LLO). These data are matched filtered through 
banks of templates that model the expected signal from a bi- 
nary coalescence of two compact massive objects with masses 
mi and m-x, resulting in triggers that pass a pre-set SNR 
threshold. We search for coincident triggers in time and tem- 
plate masses, between two or three detectors. We subject these 
coincident triggers to several tests to su ppr ess noise fluctua- 
tions (including the \ 2 test described in 02111 "). and rank-order 
the remaining coincident triggers according to their inconsis- 
tency with the background. 

We estimate the background from accidental coincidences 
by looking at time-shifted coincident triggers, as detailed in 
Section HTDI Coincident triggers that are not consistent with 
the estimated background are followed up with many addi- 
tional consistency checks, designed to identify strong but rare 
noise fluctuations. We estimate our sensitivity to GW signals 
through injections of simulated waveforms into the LIGO data 
stream which are analyzed identically to the data. 



A. Template Bank 

The templates used for this search are waveforms from non- 
spinning compact binaries calculated in the frequency domain 
using the stationary -phase approximation (SPA) l22l l23ll24ll . 
The waveforms are calculated to Newtonian order in ampli- 
tude and second PN order in phase, and extend until the 
Schwarzschild innermost stable circular orbit (ISCO). The 
templates for this single search cover a larger binary mass re- 
gion than in previous searches 111 111 , with a total mass (M) of 
2Mq < M < 35Mq and a minimum component mass of 
1 Mq. The templates are placed with a hexagonal spacing 
j25tl such that we lose less than 3% of the SNR due to using 
a discrete template bank to cover the continuous parameter 
space spanned by the two component masses. 



B. Analyzed and Vetoed Times 

The pipeline is applied to data from the first year of the 
LIGO S5 run for which more than one detector was in ob- 
servation mode. This comprises 0.419 yr of triple coincident 
data (H1H2L1), 0.232 yr of H1H2 coincident data, 0.037 yr of 
H1L1 coincident data, and 0.047 yr of H2L1 coincident data. 
In determining our upper limits, we exclude approximately 
9.5% of the data that were used to tune the pipeline lfl8ll (the 
playground data). We also exclude all the data when only the 
HI and H2 detectors were in observation mode, because of 
the difficulty in determining the background from coincident 
noise triggers in these collocated detectors (Section uTDI ). We 
make use of the (rather large amount of) additional informa- 
tion on the state of the detectors and the physical environment 
to define data quality (DQ) criteria (Appendix lAl. We use 
these DQ criteria to veto triggers in times when an individual 
detector was in observation mode, but we have reason to be- 
lieve the data were contaminated by instrumental or environ- 
mental problems. We define four categories of vetoes from 
these DQ criteria, based on severity of the data quality issue, 
and how well we understand its origin, explained in Appendix 
lAl We follow up detection candidates after successively ap- 
plying each veto category (Appendix IB1. We exclude from 
the upper limit calculation times flagged with DQ vetoes in 
the first three categories, along with triggers recorded in those 
times. This results in non-vetoed, non-playground observa- 
tion times of 0.336 yr for H1H2L1, 0.020 yr for H1L1, and 
0.041 yrforH2Ll. 



C. Coincidence Test and Clustering 

Our analysis applies a more sophisticated coincidence test 
than the one used in the past. Previously, in order for trig- 
gers from different interferometers to be considered coinci- 
dent triggers, they needed to pass a series of independent win- 
dows in time, chirp mass (M c = rf/ 5 M), and symmetric 
mass ratio (rj = m\m-2.j M 2 ). These windows were defined 
independently of the parameters of the triggers (e.g. M c , rj). 

We have replaced this coincidence test with the one ll26Tl 
that is based on the metric used in constructing a template 
bank J^^IHlllgt]. The metric contains terms necessary for 
measuring distances and determining coincidence in masses 
and time as well as the correlations between the parameters 
expected for real signal events in the 3-dimensional parame- 
ter space. This provides improved separation between signals 
and background from accidental coincidence of noise triggers, 
compared to the above independent windows. 

We have also changed the algorithm used to cluster single- 
detector triggers in our pipeline. Previously, the triggers were 
clustered by retaining the trigger with the largest SNR from 
all the templates over a fixed window of time. At present, 
we use a new method lT30tl to cluster triggers, analogous to 
the coincidence algorithm, again retaining the trigger with the 
largest SNR from a particular cluster. 
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D. Background Estimation 

As in the previous searches, we estimate the background 
due to accidental coincidences of noise triggers by repeating 
the analysis with the triggers from different detectors shifted 
in time relative to each other, forming 100 experimental trials 
with no true signals. We refer to these as time-shifted coin- 
cident triggers, as opposed to the in-time coincident triggers 
obtained without the use of time-shifts. 

This procedure is known to underestimate the rate of ac- 
cidental coincidences of noise triggers from the HI and 
H2 detectors, since they are collocated and exhibit time- 
correlated noise excursions. We therefore exclude H1H2 
double-coincident data from the upper limit calculation. We 
examine only the very strongest H1H2 double-coincident de- 
tection candidates (including H1H2 coincidences that did not 
appear in LI), and subject them to very stringent scrutiny. 
There were no H1H2 candidates that survived these checks. 
(See Section UTTIfor details.) 



E. Detection Statistic 

In this search, we employ a new detection statistic which al- 
lows us to search over a large region of parameter space with- 
out being limited by a high background false alarm rate (FAR) 
from a smaller subregion. In Ref. ifllh . coincident triggers 
were ranked by combined effective SNR (AppendixO. Here, 
instead, we use a statistic derived from the background FAR, 
as detailed in Appendix|D] The time-shifted triggers provide 
an estimate of the FAR for each in-time coincident trigger. By 
counting the number of time-shifted triggers with an effective 
SNR greater than or equal to the in-time coincident triggers' 
effective SNR, and dividing by the total amount of time we 
searched for time-shifted triggers, we calculate the FAR for 
each in-time coincident trigger. This procedure is done sepa- 
rately for different regions of parameter space with the result 
that the FAR as a function of effective SNR varies over the 
parameter space. In-time coincident triggers with the largest 
inverse false alarm rate (IFAR) are our best detection candi- 
dates. 



III. DETECTION CANDIDATES 

At the end of our pipeline we are left with a set of coin- 
cident triggers that are potential detection candidates. The 
cumulative distribution of events above a threshold IFAR is 
shown in Fig.Q] This figure shows that the loudest candidates 
in all three sets were consistent with the estimated background 
and thus were likely accidental coincidences. Thus, the search 
yielded no detection candidates, and we report an upper limit 
in Sections [Vl and [VTl 

As an exercise to prepare for future detections, we carry the 
loudest several events (such as the three loudest events that 
appear in each of the histograms in Fig. [T} through a detection 
checklist described in AppendixlBl The methods employed in 
this checklist are tested against simulated GW signals and the 



time-shifted coincidence triggers used to estimate the back- 
ground. 

Even though we know our background is underestimated 
for H1H2 coincident triggers, we reviewed the two loud- 
est H1H2 candidates using the detection checklist. In both 
of those cases, the waveforms from the two interferometers 
failed to match each other in detail, thus ruling them out as 
gravitational wave events. 

During the analysis, and prior to unblinding the non- 
playground data, an error was found in the coincidence al- 
gorithm (AppendixlBl. This caused the coincidence require- 
ment to be tighter than initially intended. It had a negligible 
effect for low mass templates, but became more significant 
at higher masses. However, since the coincidence threshold 
was selected based upon the examination of simulated signals, 
we decided to use this search in generating the upper limits 
presented here. We verified the detection candidates by re- 
running the search after correcting the coincidence test. The 
results of the corrected search did not provide any plausible 
gravitational-wave signals. 



IV. DETECTION EFFICIENCY 

We evaluate our efficiency for detecting GW signals from 
CBC during the first year of S5 as a function of mass and 
of distance to the source. This is done by coherently inject- 
ing a large number of simulated signals, called software in- 
jections, into the detector data streams. Those data are then 
analyzed with a pipeline identical to that used to search for 
detection candidates. The distribution of masses, distances, 
sky locations, orientations, and component spins is described 
in Appendix [E] The procedure for calculating the detection 
efficiency is described in Appendix[F] where a software injec- 
tion is considered to be detected if its IFAR exceeds that of 
the coincident in-time trigger with the highest IFAR. We find 
that our detection efficiencies are consistent with expectations 
from the detectors' noise spectra during S5. 

As noted in Sections U and UU we are using non-spinning 
templates to look for GW from CBC, whereas true GW sig- 
nals from CBC will have some amount of spin associated with 
the objects. Therefore, in the next two sections, we evaluate 
our detection efficiency using injections of both non-spinning 
and spinning simulated signals. 

Appendix U describes a comparison of the pipeline de- 
scribed above with one using phenomenological waveforms 
lfH [M [l5l [Tel IH. The present pipeline admits the use of 
the x 2 test 12111 . which reduces the false alarm rate at a given 
SNR threshold. Because of this, when we reduce the SNR 
threshold of the present pipeline to find the value that gives 
the same false alarm rate as the phenomenological pipeline, 
we effectively compensate for the lost signal power associ- 
ated with using non-spinning templates to search for spinning 
systems. 
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FIG. 1: The cumulative distribution of events above a threshold IFAR, for in-time coincident events, shown as blue triangles, from all 
coincidence categories for the observation times H1H2L1, H1L1, and H2L1 respectively. The expected background (by definition) is shown as 
a dashed black line. The 100 experimental trials that make up our background are also plotted individually as the solid grey lines. The shaded 
region denotes the N 1 ^ 2 errors. 



Coincidence Time 


H1H2L1 


H1L1 


H2L1 


System 


BNS 


BBH 


BHNS 


Observation Time (yr) 


0.336 


0.020 


0.041 


Component Masses 


1.35/1.35 


5.0/5.0 


5.0/1.35 


Cumulative Luminosity (Lio) 


~ 250 


~ 230 


~ 120 


(M ) 


Calibration Error 


21% 


3.9% 


16% 


Aiorizon (Mpc) 


~ 30 


~ 80 


~ 50 


Monte Carlo Error 


5.4% 


16% 


13% 


Cumulative Luminosity 


250 


4900 


990 


Waveform Error 


26% 


11% 


20% 


(Lio) 


Galaxy Distance Error 


14% 


13% 


6.1% 


A [Eq.|G2] 


0.30 


0.59 


0.45 


Galaxy Magnitude Error 


17% 


17% 


16% 


Marginalized Upper 
Limit (vr^L^, 1 ) 


3.9 x 10~ 2 


2.5 x 10 -3 


1.1 x 10~ 2 


A [Eq.|G2J 


0.30 


0.41 


0.72 



TABLE I: Detailed Results of the BNS Upper Limit Calculation 
Summary of the search for BNS systems. The observation time is 
reported after category 3 vetoes. The cumulative luminosity is the 
luminosity to which the search is sensitive above the loudest event 
for each coincidence time and is rounded to two significant figures. 
The errors in this table are listed as logarithmic errors in the 
luminosity multiplier based on the cited sources of error. 



TABLE II: Overview of Results of the Upper Limit Calculations 
Summary of the search for BNS, BBH, and BHNS systems. The 
horizon distance is the distance at which an optimally oriented and 
optimally located source with the appropriate mass would produce 
an trigger with an SNR of 8 in the 4 km detectors and averaged over 
the search. The cumulative luminosity from H1H2L1 time is 
rounded to two significant figures. 



V. UPPER LIMITS NEGLECTING SPIN 

In the absence of detection, we set upper limits on the rate 
of CBC per unit Lio, for several canonical binary systems and 
as a function of mass of the compact binary system. 

For each mass range of interest, we calculate the rate upper 
limit at 90% confidence level (CL) using the loudest event for- 
malism HHIHl], described in Appendix |Gl In the limit where 
the loudest event is consistent with the background, the upper 
limit we obtain tends toward lZgo% ~ 2.303/(TCl), where T 
is the total observation time (in years) and Cl is the cumula- 
tive luminosity (in Lio) to which this search is sensitive above 
its loudest event. We derive a Bayesian posterior distribution 
for the rate as described in Ref. 13311 . 

In order to evaluate the cumulative luminosity we multiply 
the detection efficiency, as a function of mass and distance, 
by the luminosity calculated from a galaxy population 0]] for 
the nearby universe. The cumulative luminosity is then this 
product integrated over distance. The cumulative luminosity 
for this search can be found in Table HI1 



We apply the above upper limit calculation to three canon- 
ical binary masses as well as calculating the upper limit as 
a function of mass. Our three canonical binary masses are 
BNS (mi = m 2 = (1.35 ± 0.04) M Q ), BBH (mi = 
m 2 = (5 ± 1) Mq), and BHNS (m x = (5 ± 1) M Q , m 2 = 
(1.35 ± 0.04) Mq). We use Gaussian distributions in compo- 
nent mass centered on these masses, with standard deviations 
given above following the ± symbols. 

We combine the results of this search from the three differ- 
ent observation times in a Bayesian manner, described in Ap- 
pendix iGl and the results from previous science runs |[Tll[l3ll 
are incorporated in a similar way. 

Assuming that spin is not important in these systems, we 
calculate upper limits on the rate of binary coalescences using 
our injection families that neglect spin (Appendix [E]). There 
are a number of uncertainties which affect the upper limit, in- 
cluding systematic errors associated with detector calibration, 
simulation waveforms, Monte Carlo statistics, and galaxy cat- 
alog distances and magnitudes lfl9Tl . We marginalize over 
these, as described in Appendix IH1 and obtain upper limits on 
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the rate of binary coalescences of 

^90%,bns = 3.9 x 1CT 2 yr^LicT 1 (1) 
^90%,bbh = 2.5 x 1(T 3 yr-^m- 1 (2) 
ft 9 o%,BHNS = 1.1 x 1(T 2 yr-^xo- 1 (3) 

We also calculate upper limits for two additional cases: as 
a function of the total mass of the binary, with a uniform dis- 
tribution in the mass ratio q = mVjml, and as a function of 
the mass of the black hole in a BHNS system, holding fixed 
the mass of the neutron star at ttins = 1-35 M Q (Pig. O. 



We set upper limits on the rate of these types of events that 
are two orders of magnitude smaller than the previous obser- 
vational upper limits ITTl [l3tl . although they are still several 
orders of magnitude above the range of astrophysical esti- 
mates JH 0, |a HU. In the coming years, LIGO and other 
ground-based detectors will undergo significant upgrades. We 
expect to be able to significantly improve our sensitivity to 
gravitational waves from compact binary coalescences and are 
preparing for the first detections and studies. 



Acknowledgments 



VI. UPPER LIMITS INCLUDING SPIN 

Above we have reported upper limits on the rate of merg- 
ers for different classes of objects using injection waveforms 
generated assuming non-spinning objects. We can also eval- 
uate the upper limits using injection waveforms that take into 
account the effects of spinning bodies. 

Since the maximum possible rotational angular momentum 
S for a black hole of mass m is Gm 2 /c, it is useful to describe 
the spin of a compact object in terms of the dimensionless spin 
parameter a = (cS) / (Gm 2 ). The distribution of black hole 
spin magnitudes within the range < a < 1, as well as their 
orientations relative to binary orbits, is not well constrained 
by observations. To illustrate the possible effects of BH spins 
on our sensitivity to BBH and BHNS signals, we provide an 
example calculation using a set of injections of signals sim- 
ulating systems whose component objects have a uniformly 
distributed between and 1 (Appendix|E]i. On the other hand, 
assuming a canonical mass and uniform density, astrophysical 
observations of neutron stars show typical angular momenta 
corresponding to d « 1 i34}l . In addition, the spin effects 
are found to be weak for the frequency range of interest for 
LIGO [35] so the BNS upper limits in Section [V] are valid 
even though we have ignored the effects of spin. 

Using the above injections, we obtain marginalized upper 
limits on the rate of binary coalescences of 

^ 9 o%,BBH = 3.2 x 10" 3 yr^LuT 1 (4) 
ft 9 o%,BHNS = 1-4 x 10" 2 yr^LicT 1 (5) 

VII. CONCLUSIONS 

We have searched for gravitational waves from coalesc- 
ing compact binary systems with total mass ranging from 2 
to 35 Mq, using data from the first year of the fifth science 
run (S5) of the three Initial LIGO detectors. In doing so, we 
have investigated the efficacy of searching for BBH signals 
with 2PN SPA non-spinning templates and have found them 
to be effective even at the relatively high total mass of 35 Mq. 
Additionally, we have found the non-spinning templates can 
effectively capture spinning signals with some loss of effi- 
ciency. The result of the search was that no plausible grav- 
itational wave signals were observed above the background. 
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APPENDIX A: DATA QUALITY CRITERIA 

When analyzing data from LIGO's detectors, it is impor- 
tant to know the status of the detectors at different times. We 
define data quality flags as time intervals containing known 
artifacts introduced into the data by instrumental or environ- 
mental effects. We examine the correlation between triggers 
from an individual detector and the DQ flags, and if we find 
a correlation by comparing the rate of triggers vetoed to the 
rate we would expect from fraction of time vetoed (we call 
the dead-time), we use them as vetoes. Our understanding of 
the coupling between the effect that prompted the DQ flag and 
the resulting triggers in the pipeline is measured in part by the 
fraction of the DQ flags that are used to veto triggers (called 
the use percentage). We define four different categories of DQ 
vetoes based on the above criteria. 

We categorize DQ vetoes as category 1 vetoes when we 
know of a severe problem with the data, bringing in to ques- 
tion whether the detector was actually in observation mode. 
An example case for H2 is when loud vibrations were caused 
in the detector environment in order to test the response of 
the seismic isolation systems. We categorize DQ vetoes as 
category 2 vetoes when there is a known coupling between 
the GW channel and the auxiliary channel, the veto is corre- 
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FIG. 2: Upper limits on the binary coalescence rate per year and per Lio as a function of total mass of the binary system with a uniform 
distribution in the mass ratio (left) and as a function of the mass of a black hole in a BHNS system with a neutron star mass of 1.35 Mq (right). 
The darker area shows the excluded region after accounting for marginalization over the estimated systematic errors. The lighter area shows 
the additional region that would have been if the systematic errors had been ignored. 



lated with triggers from the individual detector, particularly 
at high SNR, and when there is a use percentage of 50% or 
greater. An example is when any of the data channels in the 
length sensing and control servo reach their digital limit. We 
categorize DQ vetoes as category 3 vetoes when the coupling 
between the auxiliary channel and the GW channel is less es- 
tablished or when the use percentage is low, but we still find 
a strong correlation between the vetoes and the triggers. An 
example is when the winds at the detector site are over 30 
Mph. We categorize DQ vetoes as category 4 vetoes when 
the coupling between the auxiliary channel and the GW chan- 
nel is not well established, when the use percentage is low, 
when the overall dead-time is several percent or greater, or 
the correlation is weak. An example is when nearby aircraft 
pass overhead. We compare all of these vetoes with the times 
of hardware injections, which measure the response of the de- 
tector to a simulated gravitational wave signal, in order to con- 
firm that the DQ vetoes are not sensitive to real signals in the 
data. 



We do not analyze data vetoed by category 1 DQ vetoes. 
We remove triggers in times defined by category 2 and 3 DQ 
vetoes from the upper limit calculation. These veto categories 
significantly reduce the SNR of outlying triggers (Pig. [3j. As 
an exercise, we follow up the loudest coincident triggers after 
each category of veto is applied, including after Category 4 
vetoes (See Appendix IB1. This allows us to investigate the 
action of the vetoes by a "case study" method. 



APPENDIX B: FOLLOW-UP PROCEDURE FOR 
COINCIDENT TRIGGERS 

As an exercise, we check our loudest coincident triggers 
with a list of tests designed to see if a statistically si gnifi cant 
trigger is believable as a detection candidate. Ref. IB7I1 de- 
scribes the tests that we perform on the trigger and the data 
surrounding it. At present, our standard tests include the fol- 
lowing: We check the integrity of the data for corruption. We 
also check the status of the detectors and the presence of any 
data quality flags in the surrounding data. We assess whether 
there could have been environmental or instrumental causes 
found in auxiliary channels at the time of the trigger. We 
check the appearance of the data at the time of the trigger 
in the form of SNR time series, the \ 2 lime series, and time- 
frequency spectrograms. 

In addition, for any statistically significant candidate that 
survives the tests listed above, we plan to do the following: 
Assess the coherence between the signals recorded by each in- 
dividual detector operating at the time of the event. Verify the 
robustness of the trigger against small changes in the pipeline 
(i.e. changes in the adjacent Fourier transform boundaries or 
changes in the calibration of the data). Check the robustness 
across pipelines by employing other search techniques to an- 
alyze the same data (i.e. CBC pipelines using different tem- 
plates or algorithms designed to search for unmodeled bursts). 
Finally, we will check for coincidence with external searches 
for gamma-ray bursts, optical transients, or neutrino events. 
(This last test is for information only, as a genuine GW event 
might or might not be accompanied by other signals.) 

As mentioned in Appendix lAl we examine the distribution 
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SNR 



FIG. 3: Histogram of triggers for the H2 detector, clustered by the trigger with maximum SNR within 10 seconds, plotted after removing 
triggers occurring during times vetoed by category 1, 2, and 3 vetoes. The tail of the SNR distribution is significantly reduced by both the 
category 2 and 3 vetoes. 



of triggers after each category of veto is applied. In case there 
is a statistically significant outlier after only Category 1 or af- 
ter Categories 1 and 2 are applied, we carry out a follow up 
exercise to see if the veto that eventually rejected the event 
was rightfully applied. There are two reasons that this could 
be important. Firstly, a very strong gravitational wave from 
within the Milky Way could cause an instrumental saturation 
of the sort that we use as a veto; this kind of problem would be 
easy to diagnose if it were to occur, since the signal would be 
strong enough for us to see in the moments leading up to the 
signal-induced saturation. Secondly, we want to guard against 
false dismissal of a candidate by other kinds of vetoes, which 
can have non-negligible deadtime associated with them. Some 
of our vetoes are associated with recognizable forms of false 
signals; we check to be sure that a vetoed loud event looks like 
that kind of false signal, and not like a genuine coalescence 
signal. In the search described in this paper, there was a sin- 
gle statistically significant outlier in the distribution of events 
after the application of veto Categories 1 and 2. The follow 
up exercise confirmed that the Category 3 test that vetoed that 
event was correctly applied. 

Before unblinding the data we discovered an error in the 
computation of the template metric. This metric is used in 
the placement of the bank and the coincidence test. The error 
caused the metric distance between templates to be overes- 
timated for the higher mass signals. This has the effect of 
causing the template placement algorithm to overcover the 
higher mass region (i.e. to produce a bank with less than the 
requested 3% loss in signal-to-noise ratio). This increased the 
computational cost of the search, but did not significantly re- 



duce the sensitivity. However, this error also affected the co- 
incidence algorithm by overestimating the distance between 
triggers for high mass signals. Since the coincidence window 
was empirically tuned on software injections and time-shifted 
coincidences the impact on the sensitivity of the search was 
not significant. Consequently, the decision was taken to un- 
blind the data using the original, sub-optimal analysis in order 
to begin studying any possible detection candidates and to use 
this result to compute the upper limit (in the absence of a de- 
tection). The decision was also taken to perform a complete 
re-analysis of the data with the corrected metric to verify the 
(non-)detection statement from the original search. The re- 
sults of the re-analysis were consistent with the original anal- 
ysis and did not produce any plausible gravitational-wave sig- 
nals. 



APPENDIX C: EFFECTIVE SNR 

For this search we employ the same definition of combined 
effective SNR as was used in the BNS searches of Ref. iflTll . 
The combined effective SNR is constructed as follows. 

The single-detector SNR is produced by matched filtering 
the data against our templates. The complex output from the 
matched filter, z, is given by 

o 

where s (/)* is the complex conjugate of the Fourier trans- 
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form of the data, h (/) is the Fourier transform of the tem- 
plate, and S n (/) is the power spectral density of the noise in 
the detector. The template normalization a is given by 

oo ^ 

2 , f h(fTh(f) 

a = 4 J s n (f) df - (C2) 

o 

The z and a are combined to give the single-detector SNR, p, 
using 



From p we define the effective SNR, p e g, as: 




(C4) 



where p is the number of bins used in the % 2 test, which is 
a measure of how much the signal in the data looks like the 
template we are searching for. In the effective SNR, we nor- 
malize the x 2 by 2p — 2 since it is the number of degrees of 
freedom of this test. 

We then combine the effective SNRs for the single-detector 
triggers that form a coincident trigger into the combined effec- 
tive SNR, p c , for that coincident trigger using: 

N 
i=l 

This definition of the combined effective SNR reduces the 
apparent significance of non-Gaussian instrumental artifacts 
since it weights the SNR by the x 2 . This effectively cuts down 
on outliers from the expected SNR distribution due to Gaus- 
sian noise. In addition, we test this definition of the combined 
effective SNR using software injections and find it does not 
significantly affecting the apparent significance of real sig- 
nals. 



APPENDIX D: FALSE ALARM RATE 

Previously flUl . we defined the loudest event for the en- 
tire parameter space based on the combined effective SNR, p c 
(Appendixlcl). Since we are searching over a larger portion of 
parameter space than before, we find that the distribution of 
combined effective SNR for time-shifted coincident triggers 
varies significantly over different portions of the parameter 
space. In general, this seems to be affected by two factors. 
We see a suppression of the combined effective SNR distri- 
butions for time-shifted coincident triggers when looking at 
triple coincident triggers compared to double coincident trig- 
gers. Also, we find smaller combined effective SNR distri- 
butions for time-shifted coincident triggers in the lower mass 
regions than in the higher mass regions. 

For this search, we have decided to divide the parameter 
space into regions with similar combined effective SNR dis- 
tributions for time-shifted coincident triggers. We separate 



the triggers into different categories, where the categories are 
defined by the mean template masses of the triggers and trig- 
gers types (triple coincident triggers found in triple coinci- 
dent time, double coincident triggers found in triple coinci- 
dent time, and double coincident triggers found in double co- 
incident time). The categories for this search are given by the 
combination of three template mass regions with divisions in 
chirp mass at M c = (3.48, 7.40) Mq with trigger types given 
by H1H2L1, H1L1, and H2L1 triggers from H1H2L1 triple 
coincident time, H1L1 triggers from H1L1 double coincident 
time, and H2L1 triggers from H2L1 double coincident time. 

Within each category, the time-shifted coincident triggers 
provide an estimate of the FAR for each in-time coincident 
trigger. When we recombine the categories from the same ob- 
servation time, the FAR of each trigger then needs to be nor- 
malized by the number of trials (i.e. the number of categories). 
This normalization bestows a FAR of 1 /T with the meaning 
that during the observation time covered by this search (T), 
there is expected to be a single coincidence trigger due to 
background with a combined effective SNR at that level. 

The IFAR is used as our detection statistic and in-time co- 
incident triggers with the largest IFAR (across all categories) 
are our best detection candidates. 



APPENDIX E: SIMULATED WAVEFORM INJECTIONS 

In order to measure the efficiency of our pipeline to re- 
covering GW signals from CBC, we inject several different 
PN families of waveforms into the data and check to see the 
fraction of signals that are recovered. The different wave- 
form families used for injections in this search include Gener- 
atePPN computed to Newtonian order in amplitude and 2PN 
order in phase using formulae from Ref. 13811 . EOB com- 
puted to Newtonian order in amplitude and 3PN order in phase 
using formulae from Refs. S HI 631 |H, PadeTl com- 
puted to Newtonian order in amplitude and 3.5PN order in 
phase using formulae from Refs. i-431 6H, and SpinTaylor 
computed to Newtonian order in amplitude and 3.5PN or- 
der in phase using formulae from Refs. IVM and based upon 
Refs. J35l [M 61 S El S3 EH H; using code from 
Ref. 113011 - Each of these families except for SpinTaylor ig- 
nores the effects of spin on the orbital evolution. 

Each of these waveform families are injected from a distri- 
bution uniform in sky location (right ascension, declination), 
uniform in the cosine of the inclination angle (t), and uniform 
in polarization azimuthal angle (ip). Each of these waveform 
families are injected from a distribution uniform in the total 
mass of the system. Each of these waveform families are also 
injected uniform in log 10 D where D is the physical distance 
from the Earth to the source in Mpc. This non-physical dis- 
tance distribution was chosen in order to test our pipeline on 
a large range of signal amplitudes. 

For the SpinTaylor waveform family, each of the compo- 
nent objects' spin magnitudes are chosen from a distribution 
uniform in the unitless spin parameter a = (cS) / (Gm 2 ), 
ranging from to 1. The component objects' spin orienta- 
tions relative to the initial orbital angular momentum are cho- 
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sen from a distribution uniform on a sphere. 



APPENDIX F: CHIRP DISTANCE 

In the adiabatic regime of binary inspiral, gravitational 
wave radiation is modeled accurately. We make use of a va- 
riety of approximation techniques S l39l l40l l43l H, 14(1 l48l 
I5II l52tl which rely, to some extent, on the slow motion of the 
compact objects which make up the binary. We can represent 
the known waveform by: 



(Fl) 



where 0o is some unknown phase. For this search the func- 
tions A{t) and 4>(t) are the Newtonian amplitude and 2PN 
phase evolution respectively, which depend on the masses and 
spins of the binary. 

The template matched filtering will identify the masses and 
coalescence time of the binary but not its physical distance D. 
The signal amplitude received by the detector depends on the 
detector response functions F + and F x , and the inclination 
angle of the source 1, which are unknown. We can only obtain 
the effective distance D c s, which appears in Eq. ( IF1I ) defined 
as 12211: 



D 



D 



off 



JfJ(1 + cos 2 1) 2 /4 + Fl (cos if 



(F2) 



The effective distance of a binary may be larger than its phys- 
ical distance. 

Since the amplitude of a gravitational wave scales as the 
chirp mass M c to the five sixths power, it is convenient to 
normalize the effective distance by this, obtaining the chirp 
distance, which is given by: 



D 



chirp 



= d, 



off 



M, 



c,BNS 

Mr. 



(F3) 



where M c .bns is the chirp mass of a canonical BNS sys- 
tem. This distance is useful in evaluating the efficiency of 
the search as a function of distance since the efficiency will 
then be approximately independent of mass. 



product over distance. We combine these with the time ana- 
lyzed to calculate the posterior on the rate for the search. This 
is, assuming a uniform prior on the rate, given by lO : 



p(m|c L) t,a) 



1 + A 



(1 + mClTA) e 



(Gl) 



where fi is the rate, Cl is the cumulative luminosity, T is the 
analyzed time, and A is a measure of the likelihood of detect- 
ing a single event with loudness parameter x versus such an 
event occurring due to the experimental background, given by 
B 



AW 



-ld€L 

Cl dx 



J_dPo 

P dx 



(G2) 



The posterior (IGU assumes a known value of Cl associated 
with the search. In reality, Cl has associated with it system- 
atic uncertainties, which we model as unknown multiplicative 
factors, each log-normally distributed about 1 with errors de- 
scribed in Appendix [Hi The widths of those distributions are 
given in Table U Marginalizing over all of those unknown 
factors, and thus over Cl, gives a marginalized posterior: 



p(fi\T,A)= Pd (C L )p(p\CL,T,A)dCL 



(G3) 



where p^ {Cl) is the combined probability distribution func- 
tion for Cl given all of those unknown factors. 

The results of several experiments (e.g. different types of S5 
observing time and previous runs such as S3 and S4) can be 
combined by taking the product of their likelihood functions; 
in the case of uniform priors, this is equivalent to taking the 
product of their posteriors, allowing us to define the rate upper 
limit \i at a confidence level a by solving: 



(G4) 



where the pi (//) are the marginalized posteriors from differ- 
ent experiments calculated using a uniform prior on the rate. 



APPENDIX H: SYSTEMATIC ERROR CALCULATION 



APPENDIX G: POSTERIOR AND UPPER LIMIT 
CALCULATION 

Calculating an upper limit on a rate of coalescences in the 
loudest-event formalism requires knowledge of the cumula- 
tive luminosity to which the search is sensitive and a mea- 
sure of the likelihood that the loudest event was due to the 
observed background. The cumulative luminosity quantifies 
the potental sources of observable CBC, as measured by blue- 
light luminosity of the galaxies, which can be detected by our 
search. It is calculated by multiplying the efficiency of signal 
recovery for the search as a function of distance by the phys- 
ical luminosity as a function of distance and integrating their 



Systematic errors associated with CBC searches for GW 
signals include errors associated with detector calibrations, 
simulation waveforms, Monte Carlo statistics, and galaxy cat- 
alog distances and magnitudes. Calculating these errors in 
terms of the cumulative luminosity is described below lfl9ll . 

We refer to statistical errors associated with the efficiency 
calculation as Monte Carlo errors. Since we calculate the ef- 
ficiency as a function of distance, we calculate the error for 
a particular distance bin using the binomial formula, which 
gives an error of zero when the efficiency is zero or one, or 
when there are no injections in that bin. This error is then mul- 
tiplied by the physical luminosity as a function of distance and 
integrated over distance to get the Monte Carlo error in units 
of luminosity. 
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Calibration errors in the detectors are errors on the ampli- 
tude of the noise floor. These errors affect the amplitude, and 
in turn the distance, at which we made injections to calculate 
the efficiency of our search, since the injections we made as- 
suming a specific value of the noise floor. The one-sigma un- 
certainty in the amplitude (and thus the distance) associated 
with the calibration was 8.1% for HI, 7.2% for H2, and 6.0% 
for LI. We use these numbers to calculate the calibration er- 
ror given in Table [] in units of luminosity by combining the 
logarithmic errors in quadrature. 

Waveform errors are associated with how different the true 
signals are from what we use to measure the efficiency of our 
pipeline (i.e. the mismatch between the true signals and our 
injections). This error effectively reduces the distances in our 
efficiency calculation since we don't recover all of the power 
available in the signal due to the mismatch between the signal 
and our injections. We calculate the waveform error in units of 
luminosity assuming a waveform mismatch of 10% ll43l,l53tl . 

Galaxy errors are errors associated with our galaxy catalog 
lit] used to construct the physical luminosity. Galaxy errors 
come in two types: distance errors and magnitude errors. To 
calculate the error on the luminosity due to distance errors, 
the physical luminosity calculation is modified such that the 
galaxies' distances are increased by a factor 1 + Kj, where 
Kj is the uncertainty in the j th galaxy's distance given in the 
galaxy catalog. The galaxies' luminosities are also increased 
by a factor (1 + Kj) since the galaxies' luminosity is only 
known in terms of its magnitude and distance. To calculate 
the error on the luminosity due to magnitude errors, the phys- 
ical luminosity calculation is changed such that the galaxies' 
luminosities are increased by an amount associated with the 
magnitude errors given in the galaxy catalog. 

APPENDIX I: SPINNING SEARCH COMPARISON 

The SPA template waveforms used in this search and de- 
scribed in this paper do not take spin into account. A phe- 
nomenological template family to search for spinning black 
hole and neutron star binaries was developed in lfl2ll . referred 
to as BCVSpin, and has been used in a search of S3 data lfl3ll . 
Using both of these template banks to compute the efficiency 
of recovering signals from spinning waveforms for the differ- 
ent search methods, we find that for a comparable number of 



false alarms, SPA and BCVSpin have approximately the same 
efficiencies, implying it is not necessary to perform a search 
using BCVSpin templates in order to target spinning signals. 
The comparison of searches for spinning binaries using differ- 
ent signal models and template banks is discussed further in 
Ref. (U]. 

What is important for a search is how efficient banks are 
in picking up signals in the data. Given a large number of 
injections in the data, the efficiency is the ratio of the number 
of found injections to the total number of injections made. A 
fair comparison requires that efficiencies be evaluated for the 
same FAR. To estimate the background rate, we counted the 
number of coincident triggers in time-shifted data between HI 
and LI. 

The SNR for BCVSpin involves six degrees of freedom, 
compared to only two for SPA. As a consequence, BCVSpin 
picks up glitches more easily, and to have the same back- 
ground rate as for SPA, it needs a higher SNR threshold. (This 
problem had already been pointed out and discussed in 11211 : 
here we are seeing it in real data.) It was found that SPA with 
an SNR threshold of 5.5 and BCVSpin with an SNR threshold 
of 8 lead to comparable FARs (Fig.|4|i. 

With these SNR thresholds, we are in a position to compare 
the efficiencies of SPA and BCVSpin banks for a given FAR. 
For our purposes, an injection will be considered found if it 
had an SNR above the chosen threshold with at least one tem- 
plate in the bank, within a certain time interval around the time 
when the injection was actually made. In the case of SPA, the 
width of this interval can be chosen to be 40 ms. BCVSpin 
templates, being phenomenological, turn out to have a larger 
timing inaccuracy, and an interval of 100 ms was found to be 
more appropriate. We made 1128 injections distributed loga- 
rithmically in distance between 1 Mpc and 50 Mpc, with com- 
ponent masses randomly chosen between 1 M Q and 30 Mq 
but restricting total mass to 30 M Q , component spin magni- 
tudes 0.7 < &i < 1, % = 1,2, and arbitrary directions for 
the initial spin vectors. For the SNR thresholds of 5.5 for SPA 
and 8 for BCVSpin, in HI the efficiency of SPA came out to be 
0.93, versus 0.89 for BCVSpin; for LI the numbers are sim- 
ilar. Fig. [5] shows the efficiencies binned in distance. Hence, 
for comparable FARs, SPA and BCVSpin have approximately 
the same efficiencies showing it is not necessary to perform 
a search using BCVSpin templates in order to target spinning 
signals. 
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